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Abstract: We investigate the renormalisation group effects induced on neutrino mass and 
mixing parameters in a triplet-extended minimal supersymmetric standard model where 
a vector- like pair of hypercharge ±1 triplet superfields is added. We first rederive the 
one-loop renormalisation group equation for the effective neutrino mass operator and, for 
the case in which this operator originates solely from the decoupling of the triplets, the 
corresponding equations for neutrino masses, mixing parameters and CP-violating phases. 
We compare our results with the ones obtained previously, and quantify the importance 
of the RG induced corrections to neutrino observables by means of numerical examples. 
In the second part of the paper, we study the correlation of the model's predictions for 
the lepton flavour violating processes l{ — > £jj with the measured neutrino mass squared 
differences and mixing angles. We also emphasize the role played by the unknown reactor 
neutrino mixing angle #13 and the Dirac CP-violating phase 5. We point out that, if tan j3 is 
large, the results obtained in the commonly made approximations may deviate significantly 
from the ones following from solving numerically the relevant set of renormalisation group 
equations and using the exact one-loop formulae for the decay rates. 
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1. Introduction 

One of the most puzzling and longstanding problems in particle physics concerns the expla- 
nation of the observed fermion mass and mixing patterns. The interest around this subject 
has been renewed with the confirmation that neutrinos are massive. This motivated in- 
tense activity towards the search for an answer to the question on how neutrinos acquire 
their tiny mass. From the theoretical point of view, the idea that neutrino masses are 
suppressed by a large energy scale has become the most popular one. This is the basis of 
the well-known seesaw mechanism for neutrino mass generation. The phenomenology 
of seesaw-inspired models has been widely studied in the literature 0] with the goal of 
explaining the results provided by neutrino oscillation experiments (3), [|]. 

If the mechanism generating neutrino masses operates at high-energy scales, renormal- 
isation group (RG) effects may induce important corrections to the neutrino parameters. 
Several analyses devoted to the study of the neutrino parameter running have supported 
this expectation ||. The RG corrections to the effective neutrino mass operator depend 
crucially on the properties of the neutrino mass spectrum, on the absolute neutrino mass 
scale and on the size of the r Yukawa coupling. Therefore, the running of the neutrino 
parameters is expected to be enhanced in supersymmetric (SUSY) models with large tan/3. 
For instance, neutrino mixing may be strongly augmented by the running from high to low 
energies || so that bimaximal neutrino mixing at high scales can be made compatible with 
low-energy neutrino data by including RG corrections 0. Ultimately, with the gradual 
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increasing of neutrino data precision, even small RG effects may turn out to be important 
for neutrino mass and mixing model building. 

Within seesaw-based scenarios, the RG flow above the decoupling scale of the heavy 
seesaw mediators has to be properly accounted for, including possible threshold effects due 
to the presence of different mass scales. This has been done in Refs. || in the framework of 
the so-called type I seesaw mechanism where the heavy degrees of freedom responsible 
for the suppression of neutrino masses are singlets under the Standard Model (SM) gauge 
group. The results of such studies are model dependent since the RG effects above and 
between the thresholds depend on the Yukawa couplings of the heavy neutrino singlets 
with the leptons and Higgs fields, which cannot be reconstructed from low-energy data. 
Hence, the impact of the RG corrections strongly depends on the structure of the unknown 
fundamental couplings which is encoded in the effective neutrino mass operator. The 
impossibility of reconstructing in a model-independent way the high-scale neutrino sector 
parameters from low-energy measurements of masses and mixing angles is the main problem 
of the type I seesaw mechanism. 

An alternative version of the seesaw mechanism (usually denoted as type II or triplet 
seesaw) relies on the presence of heavy triplet states ||. Its SM version requires a single 
scalar triplet to generate mass for the three light neutrinos. In the minimal supersym- 
metric standard model (MSSM) a vector- like pair of hypercharge ±1 triplet superfields is 
demanded to ensure anomaly cancelation and holomorphicity of the superpotential. The 
analysis of the RG effects on neutrino parameters in the SM type II seesaw has been pre- 
sented in Ref. @. A more complete study covering both the SM and SUSY cases can be 
found in Ref. 0. However, the RGEs for the triplet-extended MSSM (TMSSM) derived 
in that paper differ from the ones obtained in Refs. Jl2| and ]T3| . 

The couplings of the seesaw mediators to the SM lepton doublets and/or Higgses 
may induce lepton flavour violation (LFV) in the soft SUSY-breaking Lagrangian, even if 
the mechanism which breaks SUSY is flavour blind [ [H| , |I2| ]. If large enough, such LFV 
effects can drastically increase the rates of LFV processes (like radiative charged-lepton 
decays li — > ij^f), which are otherwise suppressed to levels well beyond the sensitivity of 
future experiments. Complementarity studies between low-energy neutrino physics and 
LFV decay searches have been carried out in the context of the SUSY type I |L5| and type 
II Egl , |f6| , [ItJ seesaw mechanisms. The main difference between these two approaches is 
that the triplet seesaw is much more predictive when it comes to establishing a connection 
between low-energy neutrino physics and LFV decay searches. 

In this work we will investigate several aspects of the RG running of neutrino pa- 
rameters (masses, mixing angles and CP-violating phases) in the TMSSM, where neutrino 
masses are suppressed via the type II seesaw mechanism. The impact of the RG effects on 
predictions for the branching ratios of the LFV decays ti — > ^7 will be also discussed and 
illustrated with several examples. The layout of the paper is as follows: in Section Q we 
derive the RGE for the dimension-five effective neutrino mass operator in the framework of 
an SU(5) grand-unified model in which the heavy-triplet superfields are naturally embed- 
ded. The general form of the RGEs for the neutrino masses and mixing matrix is obtained 
in Section ||| Section |4] is devoted to the study of the neutrino parameter running in the 
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pure type II seesaw case. Namely, in Section [4,l| we obtain the RGEs for the neutrino 
masses, mixing angles and CP-violating phases (including approximate analytical expres- 
sions), pointing out the discrepancies between our results and those previously obtained 
in Ref. [pj]]. Some numerical examples are presented in Section In the second part 
of the paper we discuss predictions of the considered model for the LFV radiative decays 



j7 taking into account the latest neutrino oscillation data. We begin in Section 5.1 



by presenting the rates of these decays obtained with the help of the frequently made ap- 
proximation which neglects the RG running of the neutrino sector parameters, treats in 
a simplified way the running of the slepton mass matrices and uses a simplified formula 
for the £i — > £jj decay rates. In Section |5.2| these results are compared with the ones 
of an improved approximation in which the running of the neutrino parameters is taken 
into account and with the results obtained by solving numerically the full set of RGEs 
and by computing the decay amplitudes using the exact one-loop formulae. This allows 
us to quantify how much the approximate results of Section |5.1| deviate from the more 
accurate approaches. In particular, we show that in some cases the splitting of slepton 
masses generated by the RG running can also be important. 



2. RGE for the effective neutrino mass operator 

Let us consider a supersymmetric N=l Yang-Mills model with a superpotential of the form: 

W = — S^fc + t^Q&j + -^-$ $ 6 $ c $ d , (2.1) 

where the chiral superfields <E>j contain a complex scalar fa and a two-component fermion 
ipi, which transform as a representation Ri = R\ (g) ... <8> Rf of the gauge group G = 
G\ X ... X G n . The first two terms in the above superpotential are the ordinary Yukawa and 
mass terms and O is a non-renormalisable operator suppressed by the inverse of a large mass 
scale. Provided that higher-dimensional non-renormalisable operators only appear in the 
superpotential, then the SUSY non-renormalisation theorem still holds |l8|| . Consequently, 
the operator O can be renormalised taking into account only wave- function renormalisation. 
Using the one-loop anomalous-dimension matrices for the chiral superfields JnJ 

Y ipq = (Y^)* , (2.2) 

and following for instance Refs. [pC[| , one can write the one-loop RGE for the operator O 
as 1 : 

pabcd = Qabcf^pd + (a ^ d) + (6 ^ ^ + (c ^ d) _ (2 3) 

The quantities gu are the gauge coupling constants of the sub-groups Gk of G and Ck(i) 
denotes the corresponding quadratic Casimir invariant of the irreducible representation of 

1 From now on we will denote by X the derivative of the quantity X with respect to t = ln(A/Ao), where 
A is the renormalisation scale and Ao is a fixed but arbitrary reference mass scale. 
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We now consider an extension of the MSSM where a vector-like pair of triplet super- 
multiplets T and T transforming under the SU(3) C x SU(2) x U(l)y SM gauge group as 
T ~ (1,3,1) and T ~ (1,3,-1) is added. In a grand-unified theory (GUT) these triplet 
states may be part of the gauge group representation. For instance, in the SU(5) GUT, 
T and T are part of the 15 and 15 representations, respectively. In this case, one has 
15 = T © S © Z where S ~ (6, 1, -2/3) and Z ~ (3, 2, 1/6) under the SM gauge group. 
Below the GUT scale, the superpotential reads: 



W 

w 

W T 



W + W T + W s ,z , 

Y e e c H x L + Y d d c HiQ + Y u u c QH 2 + iiH 2 H x , 

-^(Y T LTL + XxHxTHr + \ 2 H 2 TH 2 ) + M T Ti[(ia 2 )T(icr 2 )f } , 
v2 



W s „ 



V2 



Y s d c Sd c + Y z d c ZL + M Z ZZ + M S SS . 



(2.4) 



where L (Q) is the lepton (quark) doublet supermultiplet and e c , d c and u° are the charged- 
lepton and quark singlet supermultiplets. The MSSM superpotential is denoted by Wo 
while Wt (Ws,z) contains the couplings of T, T (S, Z) with the MSSM superfields, in- 
cluding the corresponding mass terms. We adopt the SU(2) representation for the triplet 
superfields 2 
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To keep the discussion as general as possible, we will consider that an effective neutrino 



mass operator of the form [21 



W v 



■^OijtabtcdLi H 2 LjH 2 , 



(2.6) 



is also present in the superpotential. In the above equation, i,j and a, b, c, d are family and 
SU(2) indices, respectively. Notice that we neglect possible effective dimension-four Kahler 
operators of the type (Li.H 2 )(Lj.Hi)/M 2 or (Li.Hi)(Lj.Hi)/M 2 (where M is some very 
large scale) which could also give rise to an effective neutrino mass term after electroweak 
symmetry breaking (EWSB). The phenomenology of Kahler-generated neutrino masses 



(including the RG flow of neutrino parameters) has been studied in Refs. [22|. Here we 
assume that these contributions are irrelevant because they are suppressed by extra 1/M 
factors compared with the contributions of the operators included in W v . 



The RGEs for all coupling and mass parameters in (2.4) can be found in Refs. JT^] and 
1 13]. We have recomputed the full set of RGEs and found complete agreement with the 
results obtained in the latter reference 3 . Using Eq. ( |2.3[) we obtain the one-loop RGE for 



O v in (|J) as 



2 7 g 2 )if2 <T 



2 The representations of S and Z can be found in Ref. Jl^ ]. 

3 The RGEs obtained in Ref. [|13| show minor differences with respect to those of ( 



details. 



(2.7) 

See for more 



- 4 - 



T.Z 



T.T 





Hi* 



Figure 1: One- loop (non MSSM) supergraphs relevant for wave- function renormalisation of the 



lepton (left diagram) and Higgs (right diagram) doublet superfields under the superpotential (2.4). 



The one-loop anomalous dimensions for L and H\p can be computed from Eqs. ( |2.2| ) and 
(2.4), yielding the result: 



167T f L 



z 1 z 



To 9 * + \ gl ) % 



i ' 



16vr 2 7 ^ = Tr ( y e ty e + 3Y}Y d ) + 3_^ 2 - ±£ 



16n 2 4l H2 = 3Tr(F u tyj + 3|A 2 



10' 

T01 - o5 2 



3 2 
o32 



(2.8) 
(2.9) 

M « 2 «, (2-10) 

where the underlined terms are those absent from the MSSM. These new contributions to 
the wave-function renormalisation of L and Hi 2 originate from the one-loop supergraphs 
shown in Fig. [l]. It is straightforward to check that setting Yts,z = and Ai 2 = in 



Eq. (U) one recovers the MSSM RGE for O v obtained in Refs. f23]. 
After EWSB, the effective neutrino mass matrix 



m,. 



O u v 2 sin z /3. 



(2.11) 



is generated, where v = 174 GeV and tan/3 = {H®) / {H®) = i> 2 /i>i- Its RGE can be derived 
from Eqs. (p7^)-( 2.11 ), leading to: 



where 



a v = 6 Tr(YjFj + 6|A 2 | 2 - \& - 6g 2 2 , 
5 

P u = Y^Y e + 3Y T Y T + 3Y^Y z . 



(2.12) 

(2.13) 
(2.14) 



When T, S and Z have different masses, the decoupling of these states has to be performed. 
In practice, this corresponds to switching off their interactions in the RGEs. At A = Mt a 
new contribution must be added to O v due to the decoupling of the triplet fields, in such 
a way that 



O v {M T ) -»• O v {M T ) + 



\ 2 Y T 



M T 



(2.15) 



A=M T 



all quantities being taken at A = Mt- Below this scale, the running of the effective neutrino 
mass matrix follows the RGE ( |2.12 ) with Yp = and Ai = A2 = 0. If S and/or Z is lighter 
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than T, then it should be subsequently decoupled without adding any contribution to O v 
since this operator arises solely from integrating out the SU(2) triplets 4 . 

3. Running of the neutrino masses and mixing matrix 

In this section we derive the RGEs for the neutrino masses, the neutrino mixing matrix 
and the charged-lepton Yukawa couplings in the framework of the model presented in the 
previous section. We consider that at each scale A the following relations hold 

XJl m v U v = diag(mi, m 2 , m 3 ) , U\ Y^Y e U e = diag(y 2 , y 2 , y 2 ) = d\ , 

4 Y}Y T U t = diag(yf, yl yl) = d\ , (3.1) 

where U e and U v are 3x3 complex unitary matrices and rrii denotes the effective neutrino 
masses. The neutrino mixing matrix U at A is then given by: 

U = U\U V . (3.2) 

To obtain the RGE for U, we adopt the procedure first introduced in Ref. [24] for the 
renormalisation of the Cabbibo-Kobayashi-Maskawa mixing matrix, and later used for the 



neutrino case [25, £6], 27]. We start by considering the ansatz 

U e = U e Q , U v = U U R, (3.3) 

where Q and R are anti-Hermitian matrices due to unitarity of U e and U u . In order to 
determine Q, one needs the RGE for the charged-lepton Yukawa matrix Y e : 

167r 2 Y e = a e Y e + Y e P e . (3.4) 

The expressions for a e and P e are obtained by considering the one-loop anomalous dimen- 
sions Ji^ 3 and ^h[ Hi gi yen m Eqs. ([2.8D and ( |2.9[ ), respectively, and the same for the 
charged-lepton singlets 7 c 3 = 2{Y*Y^)ij — 6gf/5Sj. This yields 

i 

a e = Tr(Y e ty e + 3Y^Y d ) + 3 |Ai| 2 - ^gj - 3g 2 2 , P e = 3Y e ^Y e + 3 y}y t + 3 y\y z , (3.5) 

which, together with the diagonalisation ofYeY e given in ( |3.1| ), leads to the following result 
for Q 

W7T 2 Q ij = (P:) ij ^^f (Mi)- (3-6) 
In this equation, P' e is given by 

Pi = U\P e U e = 3(d 2 e + U\U T 4 U ] T U e + U\ Y*Y Z U e ) . (3.7) 



4 In this work we neglect the one-loop threshold corrections to the effective neutrino mass operator arising 
from the decoupling of the heavy triplet states. 
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Since Pj is Hermitian and a e is real, one has Qu = 0. Following the same procedure as for 
Q, but using now Eqs. (2.12)-(2.14) and fl3.1|), one can derive the expression for the matrix 



R: 

16vr 2 ^. = (Kh + (Pu% (i * j) , (3-8) 

where P' v is now defined as: 

Pi = Ul P v U u = uU 2 e U + 3(Ul U T 4 U V U T + Ul Y\Y Z U v ) . (3.9) 

Similarly as for Q, Ra = since Pi is Hermitian and a v is real. Finally, from Eqs. ( |3.2j ) 
and ( |3.3[ ) (and taking into account that Q is anti- Hermitian) one obtains: 

tf = -QU + UR , (3.10) 



with Q and i? given as in Eqs. ( |3.6| ) and fl3.8| ). The first and second terms on the right-hand 
side of the above equation account for the contribution to the running of U coming from 
U e and U u , respectively. 



As for the neutrino masses mj, their RGEs can be derived using Eqs. (2.12) and (|3.1j), 
leading to: 

16TT 2 mi = [a v + 2 (I*)*] m t . (3.11) 
For the charged-lepton Yukawa couplings y ei (ej = e, //, r) the RGE reads 

W7T 2 y ei = [a e + 2(P^ ll ]y e% . (3.12) 

Notice that the presence of the new couplings Ai and Yt z m a e an d P' e may affect the 
running of y e ^ T . 

4. The pure type II seesaw case 

In general, the flavour structure and magnitude of the couplings Ys,z depend on the specific 
details of the SU(5) model considered. Although not directly related with m u at tree-level, 
these couplings affect the renormalisation of Yp. It is worth mentioning that, even if one 
imposes Ys t z = Yt at e.g. A = Mq (where from now on Mq denotes the grand- unification 
scale), the RG running will deviate Y$ z from the Yt trajectory. For simplicity, we will 
restrict ourselves to the case in which the couplings I5 z are negligible when compared with 



Yt. This could for instance result from SU(5) breaking effects, as discussed in Ref. |12|. 
Therefore, from here onwards we set Yg z = 0. Moreover, in the rest of this work we as- 
sume that the only contribution to the effective neutrino mass operator O v arises from the 
decoupling of the heavy triplet states T and T. In other words, we do not address alter- 
native scenarios where extra contributions to the effective neutrino operator are present. 
Therefore, we have: 

\ 2 Y T 



(4-1) 

A=M T 



at A = M T . 
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Although O v = for A > Mt, one can define a would-be effective neutrino mass 
operator O u = X 2 Yt/M at any scale A. As a consequence of the SUSY non-renormalisation 
theorem, the RGE of this object coincides with the one given in fl2,3j). The running of the 
effective neutrino mass and charged-lepton Yukawa matrices follows the same RGEs as in 
(2.12) and (3.4), respectively, with: 

P v = Y}Y e + 3 YjY T , P e = 3Y L lY e + 3 Y T Y T . (4.2) 



Since m v is now proportional to Yt, the unitary matrices U v and Ut in Eqs. ( |3.1| ) are 
identical and 

t mlm u 2 mf A 2 w| , . 

Y T Y T = —2— > V% = i V T = -T7- ■ (4.3) 



The quantity vt is the induced vacuum expectation value of the neutral-scalar component 
of T: (T°) = vt/V%- The RGE for the neutrino mixing matrix U is shown in Eq. ( |3.1C| ) 
with Q and R defined as in Eqs. ( |3.6| ) and ( |3.8[ ). Taking into account Eqs. Q3.1| ), together 
with the fact that U = ulu u and U u = Ut, the matrices P' e and P' v are now given by 

Pl = UU 2 e U + 34 , Pe = 3d 2 e + 3Ud T uK (4.4) 

In order to obtain the RGEs for the mixing angles and CP-violating phases, we adopt 
the following parameterisation for U : 

U = K v VK a , K v = diag( e ^,e^,e^) , K a = dmg{e~ ia ^ 2 ,e~ ia ^ 2 , 1) , (4.5) 

where V'e.r,^ are unphysical phases and a.\^ are CP-violating Majorana phases. The unitary 
matrix V is parameterized in the standard way 

( C12C13 S12C13 Sl3e~ i5 \ 

V = -si 2 c 23 - c 12 s 2 3Si 3 e tS C12C23 - s 12 s 23 s 13 e lS s 23 c 13 , (4.6) 
\ s 12 s 23 - c 12 c 23 s 13 e lS -ci 2 s 23 - si2C 2 3Si3e* 5 C23C13 J 

where Cij = cos#ij, Sij = sin^j and 5 is the Dirac CP-violating phase. We identify the 
mixing angles 9 12 and 9 23 as being the ones involved in the solar and atmospheric neutrino 
oscillations, respectively, while #13 denotes the so-called CHOOZ or reactor neutrino mixing 
angle. 

Although physical observables do not depend on the phases tpi, these are crucial to 
obtain the RGEs for the physical neutrino parameters (see discussion in Ref. |^7j| ). From 



(|3.10|) and fl4j) we obtain 

idiag(0 e ,^,0 T )V + V- ^ydiag(di,d 2 ,0) = -QV + VR, (4.7) 

which can be used to extract the RGEs for the mixing angles and phases. In the above 
equation, the matrices Q and R are defined by Q = K* Q and R = K a RK*. Together 
with Eqs. (KB) and (Oh, this leads to 



2 2 

lto^Qij = (Pehj Ve ; Ve 2 j (i + j) , (4.8) 

l^ 2 R l3 = 4±^| (P„) y + {Pv)ij (i * j), (4.9) 

— mf mr A — mf 

3 % j *> 
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where 



K* P' K 

<p e ip 



Pi, K/y Pj. 



3d 2 e + 3Vd T V\ 



P v = K a {P' v f K* = K 2 a V T d 2 e V*K t 



*2 



+ 3d 2 T 



(4.10) 
(4.11) 
(4.12) 



The above results show that the RGE for the mixing matrix U does not depend on the 
unphysical phases <p e ,n,T, as expected. Comparing the present model with the MSSM, it 
becomes clear that the new contribution to the running of the neutrino mixing matrix 
comes from the first term on the right-hand side of Eq. ( |4.7| ), i.e. from the effect of the 
non-trivial running of U e induced by the presence of the couplings Yt ■ 

At this point we would like to point out some discrepancies between our results and 



those obtained in Ref. |ll|] . 



1) In Ref. |llj it is claimed that the RGEs for the mixing angles and CP-violating phases 
are independent of the Majorana phases at any scale A > My. This is actually not the 
case, as can be seen from Eq. fl4.7|) . Although the first term on the right-hand side of this 
equation does not depend on 0:1,2, the second term does (through the contribution of the 
term proportional to P in R). In fact, the dependence of the RGE of U on a\p originates 
from the running of the effective neutrino mass matrix, more specifically from the term 
proportional to Y e \Y e in Eq. (12.121). On the other hand, Q is independent of the Majorana 
phases since it is defined by P e (which does not depend on 0-1,2) and the charged-lepton 
Yukawa couplings. Hence, U e does not show any direct dependence on a\ 2. 



2) Our results for the TMSSM RGEs agree with the ones of Refs. f|] and @. However, 



we find several discrepancies with those obtained in Ref. [11] where the conventions for 
the couplings entering in Wt are the same as the ones we are currently adopting. We find 
that, in order to get complete agreement between all results, the couplings Ai,2 and Yt 
must be replaced by \f2\\ t 2 and \/2Yt in all the RGEs of Ref. [11], including the one for 
the effective neutrino mass matrix. This affects the numerical pre-factors in the second 
term of the right-hand side of Eqs. ( 4.1C| )-( |4.12 ), which in our case differ by a factor of two 



from those of Ref. [11]. The same holds for all the terms proportional to (Ai^l appearing 



in Eqs. (|2.13[) and (|3.5|) . We believe that these discrepancies might be the result of an 



inconsistent definition of the Feynman rules used in Ref. [11] for the vertices involving the 
triplet states. 

4.1 RGEs for the neutrino masses, mixing angles and CP-violating phases 

We can now use the master equation ( [4.7| ) to obtain the RGEs for the mixing angles Oij 
and CP-violating phases. For each 9^ we write the corresponding RGE in the form: 



'13 



9 ij + °ij 



b>a 



-A-ij Rab 



9-- 



£Re 

b>a 



Bt-Qab , a, b =1,2,3, (4.13) 



where and 9f j contain the contributions coming from the running of m v and Y e , respec- 



tively. The coefficients Af- and Bf^ (shown in Table |l|) are determined by solving Eq. (^ 
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s 23 e iS 
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c 23 tan 6<i3 e iS 


-s 23 tan6»i 3 e i<5 


-1 



Table 1: Coefficients Af* (first three rows) and Bfj (last three rows) which enter the definition of 



the RGEs for the mixing angles Oij, given in Eq. (4.13). 



Similarly, for the full set of phases $ = (5, a\, 02, tp e , V 9 ^) fr) 



<£> l 



b>a 



5> 

6>a 



a, 6 = 1,2, 3, (4.14) 



where <F and <3? e include the terms stemming from and Y e . The coefficients yl^ and 
are given in Table [| Regarding the neutrino masses m;, the RGEs can be obtained 
replacing P' v by P v in Eq. ( [3.11 ), leading to: 



9 

16-7T m,- 



+ 2 (V^V. + ey.; 



/??,■ 



(4.15) 



As expected, rhi does not depend on the Majorana phases a\ 2- Moreover, the third term 
on the right-hand side of the above equation is independent of the neutrino mixing angles 
and CP-violating phases. In contrast, the second term (also present in the MSSM case) 
does depend on those quantities. 

The complete expressions for the RGEs of the neutrino mixing angles and phases can be 
obtained by inserting the coefficients given in Tables [l] and |2| into Eqs. ( 4.1 3| ) and ( 4.14| ), 
respectively. In general, the final result is too lengthy to be presented here. However, 



following the procedure of Ref. [27], we will expand the RGEs to the leading order in the 
(small) mixing angle #13 . Let us first concentrate on the neutrino contribution to the RGEs 
of the mixing angles, i.e. the terms 9^ in Eq. ( 4.1 3| ) . These depend on V, on the neutrino 
masses and on the charged-lepton Yukawa couplings y\.. In view of the strong hierarchy 
Ve Vfi <C y T , we will keep only the terms proportional to y^. In this limit, and in the 
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A 12 


A 13 


A 23 
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C12S12 


C13C12S23 S13C13C23 


S13C13C23 C13S12S23 
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2 cot (9i2 


2F 3 i , 21/21 

C13C23 C13S23 


2F 32 2C12V21 
C13C23 C13S23S12 


«2 


2 tan #12 


2F 3 i 2 512^22 
C13C23 C13S23C12 


2^32 , 2F 2 2 
C13C23 C13S23 


Ve 


1 


^31 S12V22 


V32 C12V2I 


C12S12 


^13^23 <-13 i 23 c -12 


t-lS^S ( 'i3623 i 12 


in 


n 


^22 




C13S23 


C13S23 


<Pt 







V32 


C13C23 


C13C23 


$ 


Rl2 


ol3 


R23 




C23V32 ^ 3 * 2 e l<5 


523^21 


1 


C13S12S23 S13C12C13 


C23C12C13 S13C13S12 


S23C23 


O'l 


2C23V32 


2^23^22 


2 


S12C13S23 


S12C13C23 


S23C23 


0' 2 


2C23V31 


2^23^21 


2 


C12C13S23 


C12C13C23 


S23C23 




C23^32 V21 


•523^22 V31 


1 










K 13 

C13S23 





- COt 6>23 







C13C23 


tan ^13 



Table 2: Coefficients A% (first six rows) and £?f (last six rows) which enter the definition of the 



RGEs for the CP- violating phases $ = (5, ai, <%2, y e > <Pn> <Pt) given in Eq. (4.14) 



zeroth order in #13, we find 
m 3 



'13 



'12 



'23 



sin(2#23)sin(20i2) 



A 

32tt 2 Am§j 

y 2 |mie iai +m 2 e iQ2 | 2 2 . . „ 
— s 23 sm(26>i2 



, r . m,2cos(o2 — 5) rm 3 cos 5 
mi cos(oi — 0) H - ; h 



32tt 2 
32tt 2 



rAm^ 

2 |m 3 + ?ri2 e 
"12" 



'"-i 2 2 I mi + m 3 e iQl | 2 



Am|i (1 — r) 



+ s 



12" 



(1 



sin(20 23 ) 



(1 



(4.16) 
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where 



Am 2 31 



Am, 



vrii 



(4.17) 



The above expressions make explicit the discrepancies between our results and those of 
Ref. From Eqs. ( |4.16| ) it is clear that the running of the neutrino mixing angles 

does depend on the Majorana phases through the renormalisation of the effective neutrino 
mass matrix m v . The above contributions to the RGEs originate from the wave- function 
renormalisation of the lepton doublets, namely from the term proportional to YeY e , present 
below and above the decoupling scale of the triplets. Therefore, Eqs. ( 4.16 ) are valid both 
in the effective and in the full theory. Not surprisingly, the results for 9^ agree with those 
obtained for the MSSM in Ref. |2~ 



Hence, the expansions given in Ref. [11] are not valid 
above the mass scale of the triplets since they do not account for the dependence of the 
RG running on the Majorana phases, which may play a crucial role in the running of the 
neutrino parameters. 

The approximate expressions for the RGEs of the Majorana phases a\ 2 read 



a o 



vt 1 


' mym 2 s\ 3 c\ 2 


4vr 2 1 


r Amij 


y 2 r 1 




4tt 2 1 





sinfai — a 2 ) + 



sin(«i — a 2 ) + 



m 3 cos (26* 23 ) 
A ™31 

m 3 cos(26>23) 



m 1 s 12 s ai + 



m 2 c 12 s a2 
(1-r) 



2 . m 2 c \ 2 s a 2 

m lSl2 s ai + (1 _ r) 



,(4.18) 
,(4.19) 



in contrast with the result 2 ~ 



obtained in Ref. [11| at zeroth order in #13. There, the 
lowest-order term in the expansion of di i2 was found to be of first order in #13, giving rise 
to the conclusion that the RG effects on «i j2 are small. The reason why the above terms 
were not obtained in [[llj] has to do with the fact that they vanish for a\ t2 = 0, which is 
the only limit in which the expressions given in that reference are valid. 

The neutrino contribution to the RGE of the Dirac CP-violating phase 8 can be ex- 
pressed in the form [27] 



Vt 



32vr 2 6 



13 



where SY-^s and 5, 



5 U) 



d (o) 



(0) are g iven b y 
7713 sin(2^i2) sin(2#23) 



Am§! 



+ 



m\m 2 s 23 
rAm 31 



sin(ci!i — a 2 ) + 



mi sin(ai — 5) + m 2 
m 3 cos(29 23 ) 



sin(a2 — 5) 



+ 



rm 3 



sin 5 



(4.20) 



(4.21) 



m\c\ 2 sin(25 



Oil) + 



Am 2 ! 
m 2 



m\s l2 smai + 



m 2 cf 2 sin u 2 
1 — r 



+ 



s\ 2 sin(25 — a 2 ) 



(4.22) 



Again, this differs from Ref. [jll]] since does not vanish (even in the limit of vanishing 



Majorana phases) and &Y_u does depend on a.\, 2 . In the limit a\ )2 = 0, our result for 5^ 



(-1) 



agrees with the one of [11]. 

Let us now consider the charged-lepton contribution to the RGEs of the mixing angles 
and CP-violating phases denoted by and in Eqs. ( 4.13 ) and ( 4.14 ). At zeroth order 
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in #i3 we obtain 



'12 



3r Am^ y e 



ylyr + y 2 e(yl-y 2 r)cos(29 23 ) 



3Aml 1 V\ 



(vl 



2 +y 2 T 



'13 — 



y 2 



32tt 2 4 y2 
3rAmt y 2 e {y 2 T 



y 2 e){y 2 -y 2 ) 

rc\ 2 ) sin(26>2 3 ) ^ 



sin(2# 



12 J 



32^4 ' 



re 



12 



sin(2# 



2.3; 



32tt 2 4 (y 2 



yi)(y 2 -y 2 ) 



sin(2#i 2 ) sin(2# 23 ) cos 5 ~ , 



(4.23) 
(4.24) 
(4.25) 



where the final results correspond to the limit y e ^ — > 0. For the Majorana phases we find 
d| 2 — 0, at zeroth order in #13, while for 5 e we have 



■ e _ 3rAm 2 31 



yl (y 2 - vl) 



32<k 2 v\ 



{yl 



yi)(y 2 



yl) 



'13 



! 1 sin(26»i 2 )sin(26'23)sin5. 



(4.26) 



Comparing our results for the charged-lepton contribution to the running of the mixing 
angles and CP-phases (in the y e ^ — > limit) with those of Ref. plj , we find a general 
agreement. The only difference, which is the consequence of the discrepancies in the RGE 
of Y e (see the paragraph preceding Section 4.1 ), is the overall factor of 3 in Eqs. ( 4.23| )-( f4.26| ) 
which replaces the factor of 3/2 of |Q]. 

From Eqs. (|423l )-(^25|) one immediately concludes that 9\ 2 13 are mainly controlled by 
the factor rAm^/^ = y\ — y 2 , with #f 3 further suppressed by y 2 /y 2 <^ 1. On the other 
hand, #^3 i s essentially governed by Am^/vj, = y 2 — y\. Therefore, since r is small we 
expect larger RG effects on #23 than on the remaining mixing angles. This is in contrast 
with what happens for the contributions #£) shown in (4.16), where the running is typically 
enhanced for 6*12 with respect to #13,23 due to an 1/r factor present in 9\ 2 - Since the overall 
signs of #i 2j 23 an d #12,23 are the same, both the neutrino and charged-lepton contributions 
to #12,23 tend to affect the RG flow of these mixing angles in the same way. As for the 
running of #13, the main contribution to #13 comes from #| 3 which, for a given value of rrii, 
may be positive or negative depending on the values of ai,2 and 8. 

The RGEs for the neutrino masses rrii and the charged-lepton Yukawa couplings y ei 



can be obtained from Eqs. ( |3.11 ) and (3.12) with a u , a e and P' v e given as in ( |2.14 ), 



and (|4.4|), respectively. Here we focus on the running of the parameter r defined in ( |4.17j ). 
The value of r is crucial in model building since, although not affected by overall factors 
in the effective neutrino mass matrix, it is sensitive to the flavour structure of m v . At low 



energies |r(mz)| — 0.03 (see Section |4^) . 

We consider two types of neutrino mass spectra, 

Normally — ordered (NO) : m\ < m 2 < WI3 , 
Inversely — ordered (10) : 777.3 < m\ < m 2 , 



(4.27) 
(4.28) 



in such a way that r is positive (negative) for the NO (10) case. Using Eqs. ( p. 14 ), (3.11) 
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and ( |4.4| ), we can write the RGE for r as 



5. 



NO : 4vr 2 r = -rAP> 2 + (AP^ - rAP^) , (4.29) 



Am 31 



m 



10 : 4vr 2 r = AP^r - 1) + — -A- (AP^ - rAi&) , (4.30) 

with AP^ = (P^)m — (P' u )jj. To better distinguish the two main sources of RG effects, we 
express r in the form: 

r = r e + r T , (4.31) 

where r e and r T denote the terms depending on the charged-lepton and Yt Yukawa cou- 
plings y ei and yi, respectively. The contributions r e originate from the first term of P' v (see 
Eq. (fO|)) and therefore will depend on 2/e,u,r ano ^ on the neutrino mixing parameters. At 
leading order in #i 3 (and keeping only the terms proportional to y 2 ), we obtain for both 
the NO and 10 neutrino mass spectra: 

2 ( 2 

N ° 1 K = ~ 1& { ^2- V + 3r cos(20 23 ) - 2(2 - r) cos(20 12 ) S | 3 ] + 

+ r [1 + 3 cos(20 23 ) - 2 cos(20 12 )si 3 ] } , (4.32) 

2 f 2 

10 : ^ e = " 1& { t r + 3r cos(20 23 ) - 2(2 - r) cos(20 12 ) S 2 3 ] + 

+ 4(1 -r)s 2 3 cos(2fl 12 )} . (4.33) 

It is worth noticing that in the limit of quasi-degenerate neutrinos (mi ~ m 2 ~ m 3 S> 
Am^) and/or large tan/3, the RG effects on r due to r e are enhanced. For hierarchical 
(HI) neutrino masses (NO with mi <C Am 2i ), we do not expect major effects due to the 
suppression factor of r present in the second term of ( f4.32|) . Yet, in the inverted- hierarchical 
(IH) situation (10 with m 3 <C Am^) there is an unsuppressed term in Eq. ( |4.33| ) which, 
depending on the value of tan/3, may lead to a considerable running of r. 

The remaining contribution to r (denoted by rx in Eq. ( [4.31D ) does not depend directly 
on the neutrino mixing angles since it originates from the second term of P^, which is 
diagonal. From Eqs. ( 4.2S| ) we obtain 



f =-lM yH r g- r ) = _ 3Am 3i r n -r) (4 34) 

for both the NO and 10 cases. Contrarily to the results obtained for r e , the above equation 
is exact in the sense that it does not rely on any expansion nor on any special limit of the 
couplings. An immediate conclusion that can be drawn from ( fl.34| ) is that, although r T is 
always negative, \r\ decreases (increases) when going from low to high energies, for a NO 
(10) neutrino mass spectrum. Clearly, this is only true in the limit of negligible r e . 

4.2 Numerical examples 

At low energies, the neutrino mixing angles and mass squared differences are constrained by 
solar, atmospheric, reactor (KamLAND and CHOOZ) and accelerator (K2K and MINOS) 

5 Notice that the following equations are equivalent to each other. We choose to write them differently 
for the NO and 10 cases to better identify the terms which are proportional to m? 3/ATO31. 
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Parameter 


Best-fit 


3(7 


Am§! [10- 5 eV 2 ] 


7 65+ - 23 

'• D0 -0.20 


7.05-8.34 


lAmfJ [10~ 3 eV 2 ] 


2 4n+ - 12 


2.07-2.75 


sin 2 12 




0.25-0.37 


sin 2 6*23 


u - ou -0.06 


0.36-0.67 


sin 2 6»i3 


n m+ 0016 

u - Ui -o.on 


< 0.056 


r 


0.032 


0.027-0.038 



Table 3: Best-fit values (with la errors) and 3ct allowed intervals for the neutrino oscillation 
parameters from global data including solar, atmospheric, reactor (KamLAND and CHOOZ) and 
accelerator (K2K and MINOS) experiments jf|. 



neutrino oscillation experiments. The results of a global analysis Q of the data provided 
by these experiments are summarised in Table || where the best-fit values and 3a intervals 
for 9ij, Am^ and [Am^l (as well as for |r|) are presented. 

In the following, we will show some numerical examples with the aim of quantifying 
the RG effects on the neutrino mass and mixing parameters due to the presence of the 
heavy triplets T and T. We adopt the bottom-up approach, i.e. we start at A = mz 
with the best-fit values of the low-energy neutrino parameters and evolve the full set of 
the MSSM RGEs up to A = M T - At this scale, we extract Y T according to Eq. ( fO| ) and 
run the TMSSM RGEs to the GUT scale M G ~ 2 x 10 16 GeV. The RG effects in the 
neutrino mixing matrix U and in the neutrino masses mj are governed by Eqs. ( |3.10j ) and 
(3.11), respectively. We will only consider the cases of HI and IH neutrino mass spectra, 



for which the contributions to the running coming from the neutrino sector are, in general, 
suppressed with respect to what happens in the quasi-degenerate limit. In addition, only 
the results for r and the mixing angle #23 will be shown since, as discussed in the previous 
section, the RG effects induced by Yp are less important for #12,13- 

In the lower (upper) part of the left plot in Fig. |2[ we show the values of t{Mq) 
(|r(Mc)|) as a function of the coupling 7/3 (2/2) for the HI (IH) case 6 . The triplet-mass 
Mt (always given at A = Mt) is varied between 10 9 GeV and 10 14 GeV. These limits 
correspond to the solid and dashed curves (in black for tan /3 = 5 and in blue for tan (3 = 50) , 



6 The values of j/3 and 1/2 (given at A = Mg) are extracted using Eqs. (4.3) and taking appropriate ranges 
for A2 at each value of Mt- 



1 s GeV < M T < 1 0" GeV , tanp = 5 



1 s GeV < M T < 1 GgV , tanp = 5 




0.01 0.5 1 1.5 2 2.5 3 3. 

y 3 (y 2 ) 



Figure 2: Left plot: Values of r and |r| at A = Mq for the HI (lower branch) and IH (upper 
branch) cases as a function of the largest yi coupling (1/3 for HI and 1/2 IH). The blue-filled regions 
correspond to the variation of r and \r\ in the interval Mr — 10 9 GeV (black-dashed line) to 10 14 
GeV (black-solid line) and tan/3 = 5. The blue-solid (dashed) line shows the result for tan/3 = 50 
and Mt = 10 9 (10 14 ) GeV. The horizontal pink bar denotes the low-energy 3er allowed range for |r| 
as given in Table [| while the best-fit value is indicated by the horizontal dash-dotted line. Right 
plot: the same as in the left plot but for the mixing angle #23 ■ For both plots we used 9i 3 (mz) = 



and the best-fit values for the remaining parameters 
phases are set to zero. 



'21 



and |Ar 



'31 



All CP-violating 



respectively. From Eqs. (4.31)-(4.34) we obtain for the cases under discussion: 



HI : f ~ - J^r [1 + 3 cos(20 23 ) - 2 co S (2e l2 )s 2 23 ] - |j§ r(l - r) , (4.35) 
IH : r ~ - - r)s 2 23 cos(20 12 ) + |% r(l - r) . (4.36) 



For small tan /3 and i/j < 1 we do not expect large RG running effects on r since r e and r T 
are suppressed by the small couplings y T and j/j. This is true for both the HI and IH cases, 
as confirmed by the left plot of Fig. ||. For small yi and tan/3 = 5, the value of r(Ma) is 
very close to r(mz) = 0.032 (see Table ||), indicated by the horizontal dash-dotted line. 
As the couplings yi increase, r T increases and r is mainly given by the second term on 
the right-hand side of Eqs. ( |4.35| ) and ( |4.36| ) . As expected, r decreases from low to high 
energies for the HI case since f T is negative and r{mz) > 0. In contrast, although r T < 
also in the IH limit, r(mz) < and therefore |r| increases when going from mz to Mq. 
Notice that for y 3 ^ ~ 1 the value of r(Mc) is outside the 3a low-energy allowed region 
even for the largest allowed value of Mt- In the small tan/3 limit, we obtain the following 



- 16 - 



approximate results 

HI : r(M G ) ~ r 

IH : \r(M G )\ ~ |r | 



r + (1 - r ) 



a 2-i — 1 



~ (9.1 x 10~ 3 , 2.2 x 1(T 2 ) , (4.37) 



h| + (i + h)| 



» 2-, -1 

3yf 



(0.1, 4.9 x 10" 2 ), (4.38) 



where ro = r(mz)- The numbers in parentheses correspond to the estimates for Mt = 
10 9 GeV and 10 14 GeV, respectively, taking y^i = 1. 

For the HI case, the large tan 8 results are similar to the small tan j3 ones, with the 
RG correction reaching approximately 10% for small values of j/j and tan /3 = 50 (blue 
curves on the lower part of the left plot in Fig. |2|). This correction is small due to the fact 
that the term proportional to y 2 in Eq. ( 4.35| ) is suppressed by r. Instead, a large effect is 
observed in the IH limit with large tan /3 since the first term in Eq. ( |4.36| ) is not suppressed 
by r. Therefore, in this case the contribution proportional to y 2 is important even for small 
values of y 2 . We find |r(Mc)| c± 0.2 for y 2 <C 1 and tan/3 = 50. 

Let us now turn to the RG effects to 23 . From Eqs. ftL16|) and (|4.24[) we obtain 



HI : 
IH : 



?23 



723 



y 2 T l-rcos(20 12 ) + 2c 2 2 ^ ; „, ofl N 3yf 



32vr 2 



32vr 



1 — r 

,2 



sin(20 2 3) 



32vr 2 



r c 2 i2) 



3y 2 



sin(20 23 ) + ^(l 



re 



12) i 



(4.39) 
(4.40) 



which show that #23 decreases (increases) from low to high energies for the HI (IH) neutrino 
mass spectrum. This can also be seen in the right plot of Fig. [2] where #23 (Mg) is plotted 
as a function of y% (y 2 ) for the HI (IH) case. When yi > 1 the values of 02z{Mg) are outside 
the 3<t low-energy allowed interval for # 2 3 (shown in pink). Similarly as for r, the running 
effects are negligible for small tan/3 (black solid and dashed curves) and yi -C 1. The value 
of 623 at A = Mq can be roughly approximated by 



r, , x m l tan 2 8 , / M G 
9 23 (M G ) ~ fl 23 (m z ) T J- o ^ ln / 



32ttV 



32tt 2 



In 



(4.41) 



where the minus (plus) sign corresponds to the HI (IH) case. As can be seen from the right 
plot in Fig. ||, the shape of the curves is nearly the same for different values of tan /3. The 
results differ due to the second term on the right-hand side of Eq. ( |4.4l| ) , which increases (in 
absolute value) with increasing tan 8. Notice that the overall sign of that term is negative 
(positive) for the HI (IH) case, which explains the downward (upward) displacement of the 
tan/3 = 50 curves with respect to the tan/3 = 5 ones. 



5. LFV £i -»■ £j-f decays in the TMSSM 

The observation of lepton-flavour violating processes like £j — > ^7 would definitely point 
towards the existence of new sources of lepton flavour violation and/or new physics close 
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to the electroweak scale. So far, none of the aforementioned processes has been observed. 
The current upper bounds for their branching ratios (BRs) are: 



BR(/x -»■ &y) < 1.2 x KT 11 H , (5.1) 
BR(r -> e7) < 1.1 x KT 7 [§§ (9.4 x KT 8 ) , (5.2) 
BR(t -»• /i 7 ) < 4.5 x KT 8 p| (1.6 x 1(T 8 ) , (5.3) 



where the numbers in parentheses report the results for the r decays obtained through a 



combined analysis of BABAR and Belle data [31]. 

The above limits severely constrain the MSSM LFV soft SUSY-breaking mass matrices 
forcing them to be small. In the slepton sector, flavour violation can be generated in 
the presence of superpotential renormalisable interactions through RG effects. The most 
typical example is having the off-diagonal slepton masses arising due to the Dirac neutrino 
Yukawa couplings which participate in the type I seesaw mechanism of neutrino mass 
generation [14]. In spite of the huge amount of work done in the direction of establishing 



a direct connection between low-energy neutrino data and lepton-flavour violation in the 
SUSY type I seesaw, such goal cannot be achieved in a model-independent way (for a 
discussion see e.g. Ref. |33[| ). This stems from the impossibility of reconstructing the high- 
energy neutrino couplings from low-energy neutrino data. In the triplet-seesaw case the 
situation is improved since the effective neutrino mass matrix is linear on the couplings Yt 



(see Eq. 4.1 ) and therefore, in general, the flavour structure of m v is the same as the one 
of Y T . 

5.1 Approximate predictions for the ti — > £jj rates 

Let us consider the slepton soft SUSY-breaking lagrangian, 

£ soft = D m\L + e c m| c e~ ct + {H x e c A e L + H.c.) , (5.4) 

where m| and m?~ c are the SUSY-breaking masses for the slepton doublets and singlets 
respectively, and A e the trilinear terms. Starting at the scale A > Mt with universal 
boundary conditions m? = m~ c = mpll and A e = AoY e , at the scale Mt one approximately 
has pi 



(m%) i3 ~0, (5.6) 
(40v - -^MYeY^YT^ log A , (5.7) 

where i ^ j and Yt and Y e are taken at the Mt scale 7 . The above LFV terms may be large 
enough to generate observable rates for LFV processes like radiative charged-lepton decays 
h tjl- Neglecting small effects of the RG running from Mt to mz of the LFV entries 



7 In this approximation the small RG running effects on Yt and Y e between the scales A and Mt are 
neglected. 
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of (m~), the quantities on the left-hand sides of Eqs. ( p^) - (|5.7|) can be identified with 
their values at low-energies. Assuming that only the LFV coming from (m|)jj is relevant, 
and keeping the tan /? enhanced contributions to the one-loop amplitudes, one can roughly 
approximate BR(^j — > £jj) by 



48vr d a 
Gp 



CV,| 2 tan 2 /?BR(^ 



(5. 



where a and Gf are the fine-structure and Fermi constants, respectively, and BR(/i — > 
ev^Ve) = 0.1737 and BR(r — > \iv T v^ = 0.1784. The coefficients Cjj summarise the de- 
pendence of the £j —> £jj decay rate on the LFV entries of the slepton mass matrices 
and masses of the SUSY particles running in the relevant loop diagrams. In the simplest 
approximation, one has 



9l 



(m 2 



16vr 2 



, (i ^ j = e, fJ,, t) 



(5.9) 



where g<i is the SU(2) gauge-coupling constant and ms denotes a common mass scale for the 
SUSY particles participating in the process. In general, the above approximations do not 
provide an accurate result for the value of BR(^ — > £jj). For this reason, it is convenient 
to work with ratios of BRs instead of the BRs themselves. We will therefore consider the 
quantities 

_ BR(t -»• Iji) 



R 



which in the approximation (p 



" BR(/we7) 
are given by 



rt TJ — 



{m\) Tj 



BR(r 



(5.10) 



(5.11) 



i.e., depend only on the relative strength of LFV in the different channels 8 . Using the 
approximate expression for (rri~)ij given in Eq. fl5.5p , the ratios of LFV entries of the 
slepton mass matrix appearing in Eq. ( |5.11 ) can be expressed in terms of the couplings Yt 
(taken e.g. at the scale Mt) as 



(5.12) 



Unlike the seesaw type I models, Yt at the Mt or A scale can in the TMSSM be uniquely 
computed if values of the neutrino parameters are specified at mz- In this section, following 
the common procedure, we will set 



(m|) Ti 


2 


(Y T Y T ) T j 






(Yj^Yrp)^ 



{Y T Y T ) r j 


2 


(mtm v ) Tj 


2 


(YrpYrpj/ie 


A=M T 




A=m z 



(5.13) 



8 As it will be discussed in the next section, this statement is only valid in the limit of quasi-degenerate 
masses for the three slepton generations. 
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neglecting the RG running of the neutrino mass matrix between mz and Mt- We will 
assess the quality of this approximation in the next sections. 

We now use Eq. ( |3,1| ) and the parameterisation of the mixing matrix V adopted in 
Q4.5D and ( |4.6| ) to express (m 2 -)^ in terms of the low-energy neutrino parameters: 

|(m|) Me | 2 oc ^3il c\ 3 [r 2 c 2 23 sm 2 (29 12 ) + a 2 sj 3 s 2 23 + a\r\s 13 c 5 sm{20 12 ) sin(20 23 )] , (5.14) 

I /\ 771 I 

|(m|) re | 2 oc ' 2 311 c\ 3 [r 2 s\ 3 sin 2 (2fli 2 ) + a 2 s 2 13 c 2 23 - a\r\s l3 c s sm(29 12 ) sin(20 23 )] , (5.15) 

4t Urp 

| (m|) T ^| 2 oc { [ 4|r|si3 c 5 sin(2# 12 ) cos(2# 23 ) + [2 6 c 2 3 - \r\ (cos(20 23 ) - 3) cos(2fl 12 )] 

xsin(20 23 )] 2 + 16r 2 c 5 s 2 3 sin(20 12 )sin(20 23 )} , (5.16) 



where vt and r have been defined in Eqs. ( |4.3| ) and ( [4. 17 ). The above expressions are valid 



for both the NO and 10 neutrino mass spectrum with a and b defined as: 

NO: a = 2(1 - |r|s 2 2 ) ~ 2 , 6 = -2 + |r| ~ -2 , (5.17) 
10: o = -2(l + |r|s? 2 ) ~ -2 , 6 = 2 + |r|~2. (5.18) 



From the results shown in ( 5.14j )-( 5.16j ), one immediately concludes that the LFV elements 



of the soft breaking masses do not depend on the Majorana phases ai >2 Rl^] . Moreover, the 
mass of the lightest neutrino (mi or m 3 depending on whether the neutrino mass spectrum 



is NO or 10) does not have any impact on the LFV terms (mV)ij, at the one-loop level [16]. 



The independence of (mtm v )ij from the Majorana phases and absolute neutrino mass scale 
does not hold once we consider the two-loop RGEs for the soft masses. In this case, terms 
of the type 

- [y^YefYrh = -0^ |> e >tkKW , (5.19) 

K 10ft Tl{Y ^ = [mlm ^ E ml ' (5 - 20) 

^ K 0f Iy}YtY}y t]i] = ^^[mimMm^ (5.21) 



will be generated. The contributions of the form ( 5.1S| ) do depend on ai j2 and those of the 



form ( |5.20 ) and ( |5.21 ) depend on the mass of the lightest neutrino mass. However, being 



a two-loop effect, all those terms are negligible when compared with the one-loop ones. 
Hence, the quantities R T j will be mainly sensitive to 0\ 3 and 5. 

The results for the ratios R TfJi and R Te obtained within the approximation described 
above and using the latest neutrino oscillation data summarised in Table ^ are shown in 
Fig. H] as functions of s\ 3 for both NO (left plots) and 10 (right plots) cases and for the 
entire possible range of S (see also Ref. Q] ) . 
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Figure 3: Allowed regions for R TfM (upper plots) and R Te (lower plots) given in Eqs. (5.1C) and 



( 5.13 ) as a function of S13 and 5, for both the NO (left plots) and the IO (right plots) neutrino 
mass spectra. In dark green (light green) we show the 3cr (best-fit) allowed regions obtained by 
varying the CP-violating phase S in the interval [0, 2tt] and using the neutrino data displayed in 
Table [| The black solid (blue dashed) [red dash-dotted] line delimits the 3cr region for S = 
(<5 = tt) [£ = tt/4]. 



For S13 — > 0, the ratios R rfJi and R re are given by (the quoted numbers are for the 
best-fit values of #12, 623 and r, given in Table |||) 



R 

Rre 



4(l T |r|cf 2 )V 

~ r 2 sin 2 (26 12 
2 



2.3 



tan 2 (9 2 3 BR(r -> e^ T P e ) = 0.18 , 



404.0 



-18.1 (IO) 
-17.7 (NO) 



(5.22) 
(5.23) 



where the minus (plus) sign in the first equality of ( 5.22j ) corresponds to the case of a 
NO (IO) neutrino mass spectrum. Taking into account the uncertainties on the neutrino 
parameters reported in Table ||, the following 3a ranges are obtained 



R TfJL = [260 (238), 696 (751)] , R T 



[0.8,4.5] x 10" 1 



(5.24) 
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where the numbers in parentheses correspond to the 10 case. 

The fact that R TIM is for S13 = larger than R Te is due to the r 2 suppression present 
in both |(m|) Me | 2 and |(m|) Te | 2 , but absent from |(m|) TM | (see Eqs. fl5.14|) - (|5.16]) ). The 
above results also show that in the limit S13 — > the ratio R TfjL is different for the NO and 
10 cases. This is due to the fact that, although |(m|r) e(Lt | is the same in both scenarios, in 
the 10 case |(7n~) riU | is slightly larger , as can be seen from Eqs. ( |5.16 ) and ( 5.17] ). 



The parameters $13 and 5 turn out to be crucial in determining the rates of the jie 
and re LFV transitions. As #13 increases, Km 2 ^)^! and |(m|) Te | are either suppressed or 
enhanced, depending on the sign of cos 5. This can be seen from Eqs. ( 5.14| ) and ( 5.15| ) 



where, for instance, cancelations among different terms are possible for specific values of #13. 
The condition (m|)ij — > automatically implies no Dirac-type CP violation in the neutrino 

sector: the Jarlskog CP invariant J is proportional to lToa[(mtrn^)i2( m t'rn L ,)i3('rni'rn u )23\ 
and since (mim^ij oc (m/~)ij — > implies J — > |34j]. However, by itself J — > 

does not imply the absence of /i — e and r — e transitions because the vanishing of (m?^)^ e 
and (m|) re for 5 = require respectively 9 : 

1 Irl cot 6>23 sin 2#i2 , 1 H tan 2 3 sin 2#x2 fr . nr \ 
s i3 = To — \ — rr2 ' s i3 = ±o — \ — rr~2 ' V 5 - 25 ) 

2 l=p|r|sf 2 2 lzp|r|sf 2 

where the upper (lower) sign corresponds to the NO (10) neutrino mass spectrum case. This 
shows clearly that simultaneous suppresion of both [i — e and r — e LFV transitions cannot 
occur: the sign of S13 in ( |5.25| ) required to suppress the former is always the opposite than 
that required to suppress the latter. Furthermore, inserting the best-fit values of Table 



in ( 5.25 ) one finds that the value of |si3 1 ~ 0.015 for which one of these two transitions is 



suppressed is far beyond the sensitivity of future reactor neutrino experiments like Daya 



Bay [36 1, Double Chooz [37| or Reno [38|. Regarding the r^i sector, we conclude that 
(m| ) TfM shows a very weak dependence on $13 and 5 since the dominant term in Eq. ( 5.16; ) 
is proportional to b 2 cf 3 sin 2 (2#23) — 4, implying 

m L ~ — ^ — ^log . (5.26) 

IV L> Ml 16vr 2 v T M T V ' 

Moreover, the t — \x transition cannot be suppressed because the limit |(?7t~) r ^| — > would 
require S13 ~ 1, which is excluded by reactor neutrino experiments. 

5.2 Large tan/3 effects on the ratios of branching ratios 



In computing the ratios R T j in Section 5.1 several simplifications were made. Firstly, the 
running of m u from mz to Mp was neglected. Secondly, the running of (m~)ij between Mq 
and Mt was treated in a simplified way and the running of (m|)jj (including its diagonal 
entries) between Mt and mz (or - the SUSY scale) was neglected. Finally, the ratios 
R T j were computed with the help of the simplified formula ( 5.1l| ). In this section we will 



analyse scenarios in which some of these simplifications lead to incorrect results. 



9 The necessary conditions for the cancelation of (mlm^)ij have also been discussed in Refs. [Q, albeit 
in a different context. 
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We first improve the approximate calculation of Section 5A by taking into account the 
2 )■■ 



running of (m~ )jj between Mq and Mt exactly. To this end we correct Eq. (5.12) replacing 



it by 



(Y}Y, 



T ± T> T 3 



T ± T / M e 



A=M T 



(mtm 



(5.27) 



A=M T 



with the right-hand side obtained by evolving m v from mz to Mt with the help of the 
RGE ( 2.12 ). Alternatively, one can use Eqs. (|1|)-(|5T|) provided we take the values of 
all the neutrino parameters in these formulae at the scale A = Mt- We will illustrate the 



effects of this improvement by comparing the ratios |(m'~) T/J /(m?) Ate | 



and |(m|) re /(m 2 ~) 



at the scale mz obtained using Eqs. (|5.12j ) and ( 5.27 ) with the ones resulting from the 
exact numerical computation of 



at the scale mz- The latter results are obtained 



as follows. We perform numerically the RG running of m v from A = mz up to the Mt 
scale, extract the couplings Yt and run them up to A = Mq using the TMSSM RGEs. 
At this scale, we impose universal boundary conditions on the SUSY-breaking terms and 
run all the couplings and masses down to low energies. We expect deviations between the 



approximation of Section 5.1 and the improved and exact approaches to increase with tan /3 



because the running of the neutrino mass matrix is stronger for larger tan /3 values. 

In Fig. |] we show |(m?) T j/(m|) /je | 2 with j = e, (i as a function of tan /3 for the values of 
the input parameters mo, fn\j2 and Aq specified in the plots. We consider two benchmark 
values of s\ 3 {mz)'- s\ 3 = (upper plots) and S13 = 0.2 (lower plots). The results obtained 
by using the full numerical procedure are shown by the black-solid lines, while the red- 
dashed curves correspond to the approximations ( |5.27 ) . The results extracted by means of 



the approximation of Section 5.1 are shown by the black dash-dotted line which, of course, 
does not change with tan (3. 

Let us first analyse the case sis(mz) = 0. The deviations between the exact and the 
approximate results increase with increasing tan j3 due to stronger RG effects. However, 
since we are considering the case of HI neutrino masses we would naively not expect such 
large effects even for large values of tan /3 because the neutrino parameters run very little in 
this case. Although this is true for #12, #23 and r, the same does not hold for #13. Starting 
with a low-energy value #13 {mz), we have at A = Mt'- 



e 13 (M T ) ^ e 13 (m z ] 



y 2 T r + ^Jr 
32vr 2 1-r 



sin(2#23)sin(26>i 2 )ln( — ) . 

\mz J 



(5.28) 



Taking Q\j,{mz) = and neglecting the RG effects on the parameters r, 612 and #23, we 
obtain the following estimate for the value of #13 at the scale Mt 



6 13 (M T ) n -6.2 x 1(T S tan 2 /Sin 



(M T 
\m z 



(5.29) 



which reasonably agrees with the exact numerical result, even for large tan /3. Although 
small, these values of #13 at Mt may have some impact on the values of (Y^Y T )ij at that 
scale. Clearly, the effect will be stronger for larger values of tan/3 and/or Mt- The fact 
that 9i 3 (Mt) is negative leads to a suppression of \{Y^Y T ) lie \ with respect to the value 
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Figure 4: Ratios 



( m |)r Al /(™|)pe| 2 (left plots) and |(m|', Te/ v"-£^ei 
of tan/3 for s^mz) = (upper plots) and si3(m^) = 0.2 (lower-plots). The black-solid lines 



?/( TO f )^e| 2 (right plots) as a function 



correspond to the exact numerical result while the red-dashed ones were obtained using Eq. (5.27) 
with Yt taken at A = Mt- The horizontal dashed-dotted line indicates the value obtained using 
Eq. ( [5. 13 ) and the low-energy best-fit values for the neutrino parameters given in Table | (all 
the CP-violating phases are set to zero). The red-dashed lines denote the values obtained using 



Eq. ( 5.27 ) with Yt taken at A ~ Mt- Filled in yellow are the regions (delimited by black-solid 
curves) (delimited by the black-solid curves) allowed for the Mt values indicated in each plot. The 
two red-dashed lines have been obtained for the two limiting values of such interval. All the plots 
have been obtained for mg = mi/2 = 300 GeV, Aq = and assuming a HI neutrino mass spectrum 
[mi — eV) . For each value of Mt the value of A2 has been chosen in such a way that for tan {5 = 50, 

br(h -> ei) ~ 1.2 x nr 11 . 



extracted using #13 = 0. This can be understood taking into account that the last term 
of Eq. ( |5.14| ) becomes negative for #13 < 0. The opposite occurs for |(l^l^) Te | since the 
last contribution in Eq. (5.15) is now positive, implying an enhancement of LFV in the re 
sector. Since |(i^l^) T „| is practically insensitive to #13, the RG effects on this quantity 
can be safely neglected. Therefore, we expect that the values of \(Yj,Y T ) Tfl /(Yj,Y 1 



2 



and \{YjY T ) Te /{Y T \Y T ) ljLe \ 2 at A = Mt are larger than the ones predicted for 6*13 = 0. 
This enhancement should be more significant in the latter case since |(y^K r ) re | grows and 
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\{Y^Y T )^ e \ is suppressed. This is shown in Fig. || (upper plots) where one can see that the 
true values of |(m?) TU / '(mf ) , if »| (left plot) and |(m~) Tfi /(m~) liR | 2 (right plot) increase with 



i|) TM /(m|) M e| 2 (left plot) and |(m|) Te /(m|)^| 
tan/3. For tan/3 = 50, the deviation of the exact result with respect to the one obtained in 
the crude approximation of Section ^Tl] (horizontal dashed-dotted lines) amounts to about 
30% for M T = 10 9 GeV (lower black-solid line) and 90% for M T = 10 14 GeV (upper black- 
solid line) in the case of | (m? ) r ^/(m|r ) Ate | 2 . For | (m? ) re /(m|) Aie | 2 , one observes deviations 
of the order of 70% and 190%, for the same values of Mt and tan/3. Notice that the 
exact results (black-solid curves) are slightly lower than the ones obtained using Eq. ( 5.27] ) 



(red-dashed curves). This deviation is due to the running of (m'-)jj between Mt and mz- 
Considering the RGE of fn\, it can be shown that (m|)jj(m^) is approximately given by: 



(m 2 L )ij(mz) 



(5.30) 



where y ei are the charged- lepton Yukawa couplings. It is clear that this effect in (m| 
can be neglected while for (m|) TiU)Te one has: 



[m 



lJTfj,,re 



(m z ) 



16vr 2 



In 



m z J 



(m|) TMjre (Mr) 



(5.31) 



l)^ e \ obtained by solving the exact RGEs are en- 



Consequently, the ratios |(m?) r j/(m 2 
hanced with respect to those obtained in the improved approximation. Obviously, this 
effect is more relevant for large tan /3. Combining Eqs. ( 5.30| ) and (5.31) we get the relation 



(m| 




2 


(m| 


\ie 


A=m z 



l-4ln(^ 
8ir 2 \vtiz 



(m 



(m?) 



l^e 



U = V, e) 



(5.32) 



A=M T 



which explains the deviation between the solid and dashed curves in Fig. f|. 

When si^(m z ) = 0.2 (lower plots in Fig. |j), the running of the neutrino parameters 
does not affect much the quantities | (m|) T j/(m|)^ e | 2 . This stems from the fact that now 
the RG correction induced on #13 is negligible when compared with the low-energy value 
of this angle, as can be seen from Eqs. ( 5.28| ) and ( |5.29D , The deviations with respect to 
the approximate results are of the order of 10% (for the largest value of M^), as confirmed 
by comparing the red-dashed and dashed-dotted lines. Once more, the difference between 
the red-dashed and solid-black lines is due to the suppression factor shown in ( |5.32 ). 

In Fig. H we confront the ratios R T ^ and R Te obtained (for the same set of parameters 



m 



/(m 2 L U 2 



as in Fig. ^ by inserting into the simplified formula ( |5.11| ) exact values of | ( 
obtained from the RG procedure with the ratios R TfJ- and R re obtained by the exact one- 
loop calculation of the individual branching fractions BR(^j — > £jj) as e.g. in Ref. p9[ . 
Naively one would expect a good agreement between the red-dashed and solid-black curves 
since the former were obtained inserting the exact numerical results for the quantities 
I (m| )jj/(m|) /je | 2 into Eq. ( 5.10j ). However, this is not the case as can be observed in all 
plots in Fig. H The observed discrepancy is due to the fact that the formulae Eqs. ( j5lj| ) 
and (|5.9[) with the same mass ms can be a good approximation only if the slepton mass 



- 25 - 



m. = 300 QeV , m,„ = 300 GeV , A „ = , s, = 




m„ = 300 GeV , m,„ = 300 GeV , A n = , s, = 



3.2 
3 
2.8 
2.6 
2.4 
2.2 
2 
1.8 
1.6 



5 10 15 20 25 30 35 40 45 50 
tanp 

m„ > 300 GeV , m,„ = 300 GeV , A. =0 , s„ = 0.2 



( m f)rJ 

I \ m l>^ 

„ _BR(r =ffi ) 
BR(/»->ct) 

d ^ (m£m„) r J 

K T) , ~ — -f — BR(r -1 IUI, 




0.22 



0.2 



0.18 



0.14 




5 10 15 20 25 30 35 40 45 50 
tanp 

m . 300 GeV , m = 300 GeV , A = , s = 0.2 



/ 2 \ |2 

K m fW 




5 10 15 20 25 30 35 40 45 50 
tanp 



5 10 15 20 25 30 35 40 45 50 
tanp 



Figure 5: Ratios i? TM (left plots) and R Te (right plots) as denned in Eq. (5.11) for 813(7712) = 
(upper plots) and S\z{mz) = 0.2 (lower plots). The black-solid lines refer to the exact numerical 
result obtained performing the numerical RG running, full calculation of SUSY spectrum and exact 
computation of the one-loop BR(^ — > fyy). The red-dashed curves indicate the result obtained by 



means of Eq. ( 5.11 ) using the exact values of | (m 2 f )jj/(m 2 f )^ e | 2 , while the dash-dotted horizontal 



line corresponds to the approximation of Section 







spectrum is not too much split: only then can the masses of the sleptons circulating in 
loop diagrams for the processes l\ — > ^7 shown in Fig. |6] be replaced by an average mass 



771,5. I n writing Eq. ( 5.10 ) we have considered that ms is the same for r — > ijj (lj = e, fx) 
and /U —> ej. This is a valid approximation if the slepton masses are nearly degenerate for 
all three generations. 

In the TMSSM with universal boundary conditions for the SUSY-breaking terms at 
the GUT scale, the soft masses for the first two generations are approximately given by: 

{m 2 f )ii{mz) ~ml + AM? , (5.33) 

where AM? is the contribution due to the running of ( 777^)22 from the GUT scale down to 
low energies, which is mainly controlled by the terms g^\M a \ 2 (where M a is the G a gaugino 
mass) present in the RGEs of the soft scalar masses. Performing the (one- loop) integration 
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Figure 6: Example of one-loop (chargino-exchange) diagrams for the decays r — > fry (left) and 
fi — > ej (right) in the mass insertion approximation. The crossed circles denote the corresponding 
LFV soft masses. 



in the limit of vanishing Yukawa couplings, we obtain: 



AM? 



1/2 



a=l 



9i{M G ) 



i gj(M T ) 



9i(M T ) 



(5.34) 



where C[ is the quadratic Casimir invariant of the representation of / under the gauge 
group G a . The ratios of the gauge couplings appearing in the above equation are given by: 



gj(M T ) 

9 4 a (M G ) 



9l(M T ) (Ma 
8tt 2 a \M T 



9j{rn z ) 
' 9t(M T ) 



8vr 2 



\m z ) 



(5.35) 



where b a = (33/5,1,-3) and b' a = b a + 7 are the /3-function coefficients in the RGEs of 
g a for the MSSM and TMSSM cases, respectively |]l2| . The sneutrino masses of the first 
two generations (neglecting the D-term contributions and generation mixing) are given by 
m? e(j ~ ml + AM? with Ca = (3/20,3/4,0). For large tan/3, the effects of the running 
due to the large r Yukawa coupling induce a nonnegligible splitting between the masses of 
v T and v e „. In the limit of small Yp couplings, one can show that 



Vt- 



8vr 2 



A 2 

^ln 



Mg 
m z 



(5.36) 



where, for simplicity, we have approximated y T between Mq and mz by a constant. The 
term proportional to y 2 . is always negative making the r-sneutrino lighter than the remain- 
ing two. In order to account for the effect of the sneutrino mass splitting on the quantities 
R T j , we correct Eq. ( |5.1l|) with a factor r/ 



Ft ~ 



{m\) Tj 



lie 



BR(t 



F^^m^ m l s ) 



F(m 



2 2 j 



(5.37) 



where F denotes the loop function of the corresponding decay amplitude and m l s the masses 
of the non-sleptonic particles running inside the loops. 
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Deviations from the limit 77 = 1 are expected to grow with increasing tan /3 because 
the mass splitting between v e ^ and v T increases 11 with tan/3. This can be observed in 
Fig. ||. The differences between the exact numerical results for R T j (black-solid curves) 
and the ones obtained by using Eq. (5.10) together with the real values of the LFV soft 
masses (red-dashed lines) increase with tan /3. Moreover, r) depends on the initial values of 
rriQ, mi/2 an d Ao which for large tan/3 leads to a variation of R T ^ and R re in the SUSY 
parameter space, even for fixed values of Mt, A2 and the neutrino parameters. Finally, it 
is worth stressing that this effect does not depend on the way through which LFV in the 
soft masses is generated. Hence, it should also be present in the SUSY version of the type-I 
seesaw. 

In order to study the behaviour of R T j in the SUSY (7774/2,777.0) parameter space, we 
present two examples in Figs. [7| and || for tan/3 = 10 and 50, respectively. Both figures 
show the variation of R T ^ (left plots) and R Te (right plots) for sis(mz) = (upper plots) 
and s\z(mz) = 0.2 (lower plots). The blue-solid lines indicate the contours corresponding 
to BR(/j —7- &y) = 10 -12 , 10 -13 and the hatched regions are excluded by the MEGA bound 
BR(/j e"y) < 1.2 x 10~ n [^|. In cyan and light grey we show the regions where the 
mass of the lightest slepton 777^ is below 100 GeV (the dashed-dotted line corresponds to 
777^ = 200 GeV) and the lightest supersymmetric particle (LSP) is charged, respectively. 
To the left of the black-dashed line the lightest chargino mass violates the LEP bound 



777 ± 

Xi 



> 104 GeV [0] 



Fig. shows that for tan /3 = 10 (moderate RG running of the neutrino sector param- 
eters and slepton mass splitting) the ratios R T ^ and R re vary in the range of 10 % at most, 
in the whole SUSY parameter space. Moreover, the approximate estimates of these ratios 
given by Eqs. ( |5.22j ) and ( |5,23 ) are in good agreement with the exact results. In contrast, 
for tan /3 = 50 (Fig. |8|) we see that the ratios of BRs R T j are considerably enhanced for 
si3(mz) = due to the RG effects on #13, as previously discussed. Also, the variation of 
Rf^i and R re in the SUSY parameter space is now more pronounced due to the larger mass- 
splitting induced in the slepton masses (see discussion above). The observed enhancement 
with respect to the case in which all sleptons are degenerate is due to the factor 77 denned 
in Eq. ( 5.37| ) which depends on tti^ and mo, as seen in Fig. ||. 



6. Summary and concluding remarks 

It is well known that RG effects may induce important corrections to neutrino masses and 
mixing. This subject has been extensively studied in the literature in the effective-theory 
framework and also in the context of the type I seesaw mechanism. In this work we have 
addressed this problem in a supersymmetric scenario where heavy triplet states are added 

10 These deviations are negligible for the ratio BR(r — } /i7)/BR(Y — > e-f) since 77 — 1, irrespective of the 
value tan /3. 

One should recall that the dependence of r\ on the initial conditions at Mq is non-trivial. Although 
we will not address the analytical treatment of this subject here, we will show some numerical examples in 
the following. The loop functions for the various LFV operators in the two-generation limit can be found 
in Ref. {fl| 
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Figure 7: Variation of i? T/i (left plots) and R Te (right plots) in the (miw, mo) parameter space 
for Sis(mz) = (upper plots) and 0.2 (lower plots) and tan/? = 10 (the values of My, A2 and 
Aq are indicated on the top of each plot and the remaining low-energy parameters are taken at 
their best- fit points). The mass of the lightest slepton is below 100 GeV inside the regions 
filled in cyan and the black dashed-dotted line corresponds to = 200 GeV. To the left of the 
black-dashed line the lightest-chargino mass is below the LEP bound. The black-hatched region is 
excluded by the MEGA bound BR(/i — > ej) < 1.2 x 10~ n 28 and the blue contours correspond 
to BR(p — > ey) = 10~ 12 , 10~ 13 . The variation of R T j follows the colour bars shown on the right of 
each panel. 



to the MSSM field content. We started by obtaining the RGE for the effective neutrino 
mass operator in an SU(5) model in which the triplet superfields can be accommodated 
in a 15-dimensional representation. The general expressions for the RGEs of the neutrino 
mixing angles, masses and CP-violating phases were also derived. Taking the pure type II 
seesaw case, we have analysed both analytically and numerically the effect of the couplings 
Yt on the RG flow of the neutrino masses and mixing. Our results can be summarised as 
follows: 

• We have pointed out some differences between the present results and those previously 
obtained in Ref. [11]. Apart from discrepancies in the RGEs, we have shown that for 
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Figure 8: The same as in Fig. |?] for tan/3 = 50. Inside the light-grey regions the LSP is charged. 

energies above Mt the RG flow of the neutrino mixing angles and CP phases does 
depend on the Majorana phases. 

• The RG-induced effects on the neutrino masses and mixing angles due to the presence 
of the heavy triplet states become more relevant as the size of the couplings Yp 
increases. These new effects are not sensitive to the value of tan /3 and may be equally 
relevant in a non-SUSY case. The running contributions controlled by Yt are more 
important for the parameter r (the ratio between the solar and atmospheric neutrino 
mass squared differences) and the atmospheric neutrino mixing angle. Regarding the 
running of the CP-violating phases 5 and a± t 2, we have shown that the RG effects 
induced by the couplings Yt are negligible. 

• Within the bottom-up approach, we have worked out some numerical examples (for 
both HI and IH neutrino mass spectra) with the purpose of quantifying the running 
effects on r and #23- We have shown that if yi ~ 0(1), then the values of #23 and r 
at the GUT scale are outside their present low-energy 3a intervals. This means that 
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type II seesaw-based models for neutrino masses and mixing with large Yp couplings 
should not predict neutrino mass and mixing parameters at a high scale in agreement 
with the low-energy data and, therefore, a consistent RG analysis is demanded. 

The second part of this work has been devoted to the analysis of the LFV charged- 
lepton radiative decays — > £j"/ and their connection to low-energy neutrino data in 
the framework of the TMSSM with universal boundary conditions for the soft SUSY- 
breaking terms. After having updated the approximate predictions for the ratios of BRs 
R T j = BR(£j —> £7-7) /BR(// —7- &y), we compared them with the exact numerical results. 
The predictions depend on the value of tan /3 and, of course, on the value of the yet unknown 
neutrino parameters #13 and 5. To summarise the results of this part, let us review the two 
extreme limits of low and high tan/3. 



Small tan/3 

If tan j3 is small, the RG-induced effects on the neutrino mass and mixing parame- 
ters can be safely neglected, and the ratios of branching ratios R T j (j = e, fj,) only 
depend on two parameters: 5 and #13. Furthermore, the approximate results are in 
good agreement with the exact numerical ones. For #13 = 0, we have shown that 
the 3<j allowed interval for R T ^ is R Tfl = [260 (238), 696 (751)] (where the numbers 
in parentheses refer to the 10 neutrino mass spectrum case). This means that, if 
BR(/i — > ej) ~ 1CP 11 (close to the present upper bound), then r — > fij could be 
observed with a future sensitivity of 10 -9 , reachable at the SuperKEKB upgrade 
and SUPERB [43] flavour factories. However, if the bound on BR(/j — > ej) is lowered 
to 10~ 12 , then an observation of r — > /ry at the level of 10~ 9 or above would be in 
conflict with the predictions of the model for small tan f3 and #13 = 0. For the r — > cy 
decay, its BR is too small to allow for its future observation. 

Considering now values of #13 close to the experimentally allowed upper bound (#13 ~ 
0.2), the predictions for R Tfl show that an observation of r — >• iJ/y above 10" 9 would 
exclude the TMSSM with universal boundary conditions for the SUSY breaking terms 
at the GUT scale. This stems from the fact that for #13 > 0.1 (and considering 
BR(^ — >• cy) < 10~ n ), BR(r — >• /ry) < 10~ 10 . For the r — > decay, its observation 
with BR(t —t- &y) > 10~ 9 would also be in conflict with the TMSSM if #13 is close 
to its upper bound. Nevertheless, for 6*13 ~ 0.015 flavour suppressions in the /ie and 
re slepton masses may occur and, under special conditions, all the three LFV decays 
could be observed. 



• Large tan/3 

We have concluded that if tan /3 is large the values of R T j may deviate considerably 
from the approximate ones. Moreover, in this tan /3 regime the splitting between 
the third generation slepton masses and the remaining two introduces a non-trivial 
dependence of R T j on the soft SUSY-breaking masses. This is in contrast with the 
low tan (3 limit in which the slepton mass splitting is much less pronounced. We have 
shown that the running of the unknown mixing angle #13 from the electroweak to 
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the heavy-triplet decoupling scale may affect the relative magnitude of various LFV 
soft slepton masses, especially in the case of a very small #13. In short, this means 
that ignoring the RG running in estimating the rates of the LFV processes on the 
basis of the neutrino data may lead to misleading predictions for the rates of the LFV 
decays. Therefore, for large tan /3, the only way to obtain reliable results is to perform 
a complete numerical calculation. In particular, R Tfl may be enhanced by a factor of 
four (or even larger) when going from tan f3 = 10 to tan f3 = 50, for the same set of 
initial conditions (see Figs. and ^). Consequently, if tan/3 is large (and #13 is close 
to zero) BR(r — > fi'y) can reach the value of 10~ 9 (and, therefore, be experimentally 
detectable) even if BR(/i — > e-y) is as low as 10 -12 . Still, the simultaneous observation 
of /j, — > ej and any of the two radiative r decays would strongly disfavour the present 
scenario for large #13, even in the high tan/3 regime. 

In conclusion, we have shown that the RG corrections in the framework of the TMSSM 
may be important when making predictions regarding neutrino masses and mixings. These 
radiative effects become more relevant when the couplings between the lepton doublets 
and the heavy triplet are large. On the other hand, provided that tan f3 is large, important 
effects may be observed on the LFV decay rates even if the couplings Yt are small. In such 
a situation, the knowledge of all low-energy neutrino parameters and SUSY mass spectrum 
is crucial for an accurate prediction of the LFV decay rates. Nevertheless, the approximate 
results for the ratios of BRs in distinct channels are reliable if tan /3 is not too large. 
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